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ABSTRACT 


XI i 


In the oresent work, the combined finite element- 
transfer matrix (^ETM) method has been developed for the vibra- 
tion analysis of cable networks and beam grillage. Firstly, 
the system is divided into elements and their transfer matrices 
are obtained from the closed form solutions of their equations 
of motion. From these transfer matrices, stiffness matrices 
are developed to obtain the transfer matrix for the repeated 
section using finite element method (FaO. This transfer 
matrix is then used to derive the transfer matrix for the whole 
system. Thereafter, the boundary conditions are applied to 
obtain the characteristic equation, roots of which give the 
natural frequencies of the system. 

For cable networks, the transverse vibrations of the 
cables are considered. Cases of the cable networks, orthogonal 
and non-ortho gonal are dealt with. Transverse vibrations of the 
beam grillage in bending and combined bending-torsional modes 
with various types of boundary conditions are analysed. 

It is conjectured that FETM utilizes lower memory and 
has a simpler formulation than FB/l for the class of repeated 
system. It has the same advantage of obtaining formulation 
from the exact solution of its elements, as is for transfer 
matrix method. 
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INTRODUCTION 

The finite element method (pll') is the most powerful 
tool hence vaoely used for the vibration analysis of structures. 
However, -^^or complex structure, the method reou? res a large 
number of nodes resulting in a very 3erge comoutor memory 
reoui recent. 

^’okainish [l] used rhe combined finite element-transfer 

matrix method (PET'') in the study of the dynamics of rectangular 

pjatos. In his method, the siize or stiffness and mass matrices, 

is ecual to the nuirber of degrees of freedom of one strip, so the 

frecuency determinant for a clamped-clamped pl<ate is 18 x 18 by 

the FETi- as compared to a 108 x 100 matrix eigen value problem 

obtained by using the standard PS/ with same num^ber of nodes. 

iVcDaniel and Eversole [?] have proposed a similar approach in 

oH-ef 

treating a stiffened plate structure. Sankai and Hoa 
an approach, in v;hich an extenoed transfer matrix relating the 
state vectors which consists of state variables (displacements 
and forces) and their derivatives with respect to frequency was 
used. Ohga and Shigematsu [4j have also used FETM for bending 
and buckling problems and have developed a technique for treating 
the structures with intermediate conditions. 

The cable networks are used for covering large column- 
free areas. This is done due to their 'light weight and high 
strength resulting in reduced cost per unit covered area. It 


also adcis to apsth<»Tac senso. The cables nave neqlic<ible fJ ex- 
'jral rigidity and sen^e onnerily as the tension merabers. 'Tonoi- 
derable "-ork has done in studying Crshle networks subjected 

to static loading whereas relativol/ a little vrork has been done 
in the study o? dynamic resoonse of cable networks i.bj more so 
due to the comolexity of the shaoe. The static and dynamic 
analyses of the cable netw'orks have been done by modelling as 
discrete systems. This apuro^cn loads to a set of simultaneous 
algebraic and o’lfferential eauations respectively and it becomes 
quite involved in terms of computer time and memory as the 
number of equilibrium equations are generally large. The memory 
requirement also bscone very large if consistent mass-matrix 
approacn is used in the dynamic analysis. 

The finite element method with lumped mass approach has 
been used by Choudhari [6]. Singh [7] has used the F3/' with the 
consistent mass-matrix approach. In his work, a cable between 
two joints has been taken as one finite element for the analysis 
of the cable networks. 

A qrillage/grid is a structure consisting of two sets 
of orthogonal beams rigidly connected at the intersections. 

Some of its applications are in building floors (concrete slabs 
embedded in grillage), ship decks, bridge floors, missile ground 
facilities, electronically steerable radar systems etc. Several 
studies [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18] have been 
reported of the static and dynamic analyses of the beam grillage. 
One of them, namely the discrete approach is suitable for the 
analysis of such systems but is laborious and time consuming. 
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Furthermore if one uses consistent mass-matrux aoproacli, it can 
be solved for higher order grillage utilizing largo raemory 
computers. 

The earlier ivork on the vibration of grillages is Icnov/n 
to bo that of 3ox ano Denko [o] \dio formulated the problem in 
the matrix form, suitable cor the use of digital computer. 
Ellington and ''cEallion [9j have taken into account the discrete 
nature of the beams, by (a) assuming the moss of beams as concen- 
trated at the intersections (thereby reducing the problem to a 
finite-degree of freedom system), and (b) neglecting the tors- 
ional stiffness of the beams. A more general solution was 
presented by Thein "/oh [10j V'/ho extended the method of Ellington 
and f’cCallion [9j by taking distributed mass of the beams, 
instead of lumping the moss at the joints. Such a study is 
necessary to assess the accuracy of the approximate solutions 
obtained by lumping the masses. However, due to complexity of 
the formulation, the solutions obtained are for bending deforma- 
tion only. Thein Wah’s [10] work was extended by Chang and 
Michelson [11] in which the complicated equations were simplified 
with the application of Rayleigh-Ritz method. 

Several studies have been made in which the grid is 
replaced by an equivalent orthotropic plate. For its static 
analysis, Timoshenko [12] considered the grid and gave the para- 
meters of an equivalent orthotropic plate in terms of the known 
parameters of the grid. Using Timoshenko's results, Tarapore- 
walla [ 13 ] found static deflection of the grillage by considering 
the plate analogy. Renton [14] gave the theoretical basis for 
the relations given by Timoshenko. He took the advantage of 



tho i'9sultin'.') regularity oT the grids, consisting of a lirqo 
nur-ber of heairs. 'le used tne coiubinobion or ’’finite differencos”, 
and oifferential calculus, the two being related by Taylor’s 
expansion. This technicue gives dif-Perential equation for the 
deflection of an eoui valent continuum. In this further work, 
lenton [ibj gave first and second order continuum approximation 
for various 'om.s of soaco grids. In the paoer of Cheng Lloj, 
the results for the frequency parameters for a 3 x 3 square 
uniform grid are given, 'vhich have been obtained from eoui valent 
orthotrooic plate [l7j. The results compare well with the results 
obtained from various dd scmte methods. Sharrna [isj has done 
the dynamic analysis of simnly supported beam grillage based on 
rational prooosition of an orthotrooic plate analogy, formulated 
on the basis of energy equivalence, thus replacing the Timoshenko’ 
plate analogy which is an intuitive proposition. Sharrna ['is] has 
further extended the proposition to the case of grillage with 
clamped edges but the results are limited to the fundamental 
frequency only. 

In the present work, FETIA method is used for the free 
vibration analysis of cable networks and orthogonal beam grillage. 
Newton-Raphson modified technique is used to determine its 
natural frequencies. 

In Chapter 2, the free vibration analysis for the cable 
networks using FETM method is given. The problem formulation is 
done for the cables lying in two directions, orthogonal and non~ 
orthogonal. Chapter 3 deals with the free vibration analysis 
of clamped orthogonal beam grillage in (a) transverse bending 
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rporio oiily one’ (b) co’^'ib: ned bonUi ncf-xorsional snodo. Method of 
solution is cisciosed in Chapter -!■. Results for various cable 
networks and bea^ gr5I}age aro presented in Chapter 4 and 
compared mth earlio'^ results, '.vherever possible. The conclu- 
sions derived fron the nresent work are oiven in Chapter 5. 


0 



0 
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The geopotrv ora cable net'v'ork with rectan: 7 Ular or 
SQuare ooundaT'y can be consideroc' to be consisting of repeated 
sections, "iiare -c^e'rceob tne cables in the nei ghbournood of 
the boundaries and hence r.'ethod is euite suitable for its 

dynamic analysis. Tq deal mth the transverse vibrations of 
the cable network ’nth ?5Tb one needs the transfer matrix (Ti') 
for a sinale cable element which in turn requires its ecuation 
of motion and its solution. In Iho succeeding section, the TK 
for a vSimqlo cable is obtained from its equation of motion. 

From the Tfi of a cable element, stiffness matrix for the cable 
element is obtained. Then elemental matrices and state vectors 
arc assembled to obtain the global stiffness matrix and state 
vector for a repeated section. From the global stiffness matrix, 
T7.' relating state vectors on right and left sides of the repe- 
ated section is obtained. Because of its repeated character- 
istic, the system Tfl is obtained by multiplying the TM^ corres- 
ponding to number of reoeated sections and finally vath TK of 
the left and right hand boundaries. . Applying the boundary 
conditions of the system the characteristic equation is arrived 
at, which is then solved for its roots. In the present section 
cables lying in two directions are considered. Directions for 
networks with cables in two directions are denoted by X and Y. 
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ri ^ , f ■'T'rr"' > \T'-'T’'’ c •> A -'A AT "'T a 


-OA’ ? cTbJe si'ovn in ?i jure 1 with masb oe : unit length 
i{X), tension T(X) an" L-’'ansvGrsG loading c(X) oor unit length, 
hi? G'-adtion of no cion is 


; [T(X; -I =r rn({) 


( l '.1 } 


O t 


v'hore 7 is cransverso disnlacenent of nho cable. 

gonsic’ar tie cable cl enent shovsni in Fi-iure 1. The TM 
for the cable is liven by eauation [iQJ 


(2.2a) 



r 1 1 

cosp p S3np 


f 

= 

1 


i 

-ps'^np cosp 



L J 





or 


where 


= [Tj 


n = '-l/l, p = Ml/C, F = f^/T, c = (T/m) 

J = length of the cable 


1/2 


(2.3) 


= vertical component of the force acting on the cable 
c = wave equation constant 
0 ) = circular frequency of harmonic motion 


and 


{Z}^ and ^ are state vectors consisting of displacement 

and force at the station i and i-1 respectively and [t] is TM. 


2.1.1 TRANSFER MATRIX FOR THE CABLE LYING IN X DIRECTION 

Let 1 , T and m be the length, tension and mass per 
unit length respectively for the cable lying in X direction. 
Then its TM relation is given by 



8 



1 



F.-1 


I 


1 




Fig. 1 A cable element 

(Non -dimensional quantities) 


n+1 




0 Node number 

1 I Element number 


Fig. 2 Finite element mesh of a repeated section. 
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or = [ 

whore r] = 


-1-i ^ .-1 


’VX^ 'C’ ^ ^1^,/c 


(l. In; 

(2.4b) 

(2.5) 


■o’ce 'ch't tile see to I'ecbo.' o~ a-iy cable lylny in 'C direction is 
identifier, by subscriot 1. 

"^.1.^ ’ \T''.'y -21 12 2'\.‘.L2 Ltl’v'S lU Y DI12CTI0\’ 

Let ]y, en<j rn^ be tne lengtn, tension and nass per 
unit length respectively for the cable lying in the Y direction. 
Tnen the PI relation is given by 




sinp^hy 


sinp^riy cosp^T)^ 


(2.oa) 


0-1 


or 


czl^ = [Tg] 


(2.5b) 


r) 2 

where r] = ’ 


r^, = 


T /T 
X y 


m 


m /m . T) =1 /I 
X Y V Y X 


(2.7) 


The state vector of any cable lying in the Y direction is iden- 
tified by subscript j . 


2.1.3 TRANSFER MATRIX FOR A REPEATED SECTION OF CABLE NETVIORK 


To derive the TM relating the left and right hand state 
vectors (displacOTient and force) of a repeated section Figure 2, 
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rst b J s .v^'] bo o' ' ^ b’io> s‘5'iffno'5s '"■Titrix ^o- tne 

coble oJo^ent. ''o ''o so, '^rr- . ;ns (C. 'o) an'"’ (O.o?) con be 
“rons -on~>oct into 


r ^ 

! 

f = 

C 1' 3 . 

- - - 

5'^ 

sTrTp 


H 

1 

lb J 


"" ’?in 

SI no 


1 

y 

J 



- 

-J 




= [g] 






(2.0a) 


(2.8b) 


'Yhero IS the stiffness 


toi\ of the cable element in X 


direction. 

Similarly foi the coble 'Element in Y direction 



cosn . I . 

-J- ...J.J. 

Tj sihe^'T)^' 
n . , 

J. . . J. 

f: c'^'rTyi^ 




b 

_ ‘JL „-J 


h 

rj sinp^T]^ 


j-1 

P-, COSp^T]^ 


T] I 

r^ slnp^T]^ 


^ J 


(2.ba) 


or CF),^ = [ly {T)}^ 


(2.9b) 


where [K^j Is the stiffness matrix for the cable element in Y 
direction . 

To develop the TM for the repeated section, shown in 
Figure 2 from those of its stiffness matrix, a numbering system 
has been adooted similar to that of F3\ . The global stiffness 
matrix is obtained from those of elements in x and ^ directions 
on the lines of FH/*. So there will be a total of 3n + 2 nodes 
and 3n + 1 elements for n cables lying in the ^ direction. 

The procedure of assembling the element matrices and 
state vectors is based on the reoui reinent of ’’compatibility” at 





'.he e.i enent node's. ti-jIs iTpana that nt the nodes i\/hero elements 
ere connected, tne vajuoCs) or tne unhno'm nodal Qecireo(s) of 
freeo’on or vnriablG(s) isCare) tne same for all che elements 
joining at the nodeo In cho oresent oroblem, the generalised 
di splace’^on t is ■*! (non-di'"'ensional displac an ent). '7hen the 
generalised displacements are 'Tiarched at a common node, the 
nodal stiffness and i^orcos oT the elements meeting at the node 
are added to obtain net stiffness and net forces at that node® 
[n view of tne above procedure, the global stiffness matrix is 
obtained as 

1 '^...n |n+1 1 n+2 n43...2n+1 I '^n+2 I 2n+3 '^n+4. . .3n+2 1 


L ‘1^ I i 

(nxn) I I 


(nxn ) 


LOj 

(nxn) 


(nxn ) 


C'aJ 

(nxn ) 


C'V 

(nxn) 


%-f2 

:^n+3 1= 


[0] 

(nxn ) 


^ ^ 


rv 

(nxn ) 


(nxn) 


_^2n+j_, 

h 

'H 

204*3 

2n+4 

• 

m 

m 

"H 

3n+2 


(2.10a) 


or [Kg] 


(2.10b) 


where [Kg] is global stiffness matrix. 
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-c nati on (f.lO) rol.’ces Lne forces and disolacoponfcs in 
tne finite elG^'’ent rosn of rhe repoaced section. The boundary 
conditions of nodes (n+1 ) and (rn+2) 

are niven by 


ri 


n4 1 




’n+? 


0 


(?.11 ) 


■^he qlobal stiffness natrix is oartitioned as indicated 
by dotted lines in onuat3 on (^.lOa). It rnay be noted here that 
all I \J matrices a-^'e symmetr] c, the tv/o corner matrices are nulJ 
matrices, and diagonal matrices. 

Equations (p.lOa) are rewritten as 

4- {r]}^ = £-?}^ (2.12) 

DApJ + LA 3 ] + t'V ^ 1)3 = CO} (2.13) 

C-A^] Cn}^ 4- D'V Ci)}3 = CF}3 (2.14) 

where {iil,, and {ti }3 are displacement vectors at left hand, 

right hand and at the joints of finite element mesh of the 
repeated section respectively; C-p}^ and Cf }3 are force vectors 
at left hand and right hand of the section under consideration 
respectively. 

To derive the TM fPr the repeated section, one needs 
to determine {'n }3 and CF }3 in terms of Ct)}^ and C F . This can 
be accomplished by eliminating {rj}^ from the equations (2.12), 
(2.13) and (2.14) and is done as fo31ov/s: 
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From equation (2.13) 

;t,) 2 = -LAjj'' rA2J - [Ajj-’ SmIj (2.15) 

Substituting this in equation (2.12) one gets 

0^3 = [T^j (ri!^ + {Fl^ (2.16) 

where [T,j = [TjjCDAiJ - 

and [Tjj = [ITAjjiAg]-'' fA^J]-'' (2.17) 

Substituting It )}2 eauations (2.15) and (2.16) in 

equation (2.14), one gets 


= [T 3 ] Cnl., + [T^] 

where [T 3 ] = ( DAgd - [ AjJ-l [A^j ) [T, ] 

-TA^JEAg]-' TAj-J 

and [T 4 ] = (iTAjJ - lA^JLAjJ-'' lA^J) [T^I 

Equations ( 2 . 16 ) and (2.18) are combined to obtain 



or iZ}^ - [T] iZ)^ 

where [T] is TM for the repeated section. 


(2.18) 


(2.19) 


(2.20a) 

(2.20b) 


2.2.1 TRANSFER MATRIX FOR ORTHOGONAL CABLE NETWORK 

Take the case of an m x n rectangular orthogonal cable 
network shown in Figure 3, where m and n are number of cables in 






15 


'( and Y directions resooc tively. The cables lyang in Y direction 
are divided into 2(n-!-1 ) eoual fields, that is, portion of each 
cable betv/een two adjacent joints comprises of two fields. In 
the case of cables lying in the Y direction, the portion of 
cable between two adjacent joints and at the boundaries is 
treated as single field. It is to be noted that equations (2,8a), 
(2.8b) and (2.9a), (2.9b) hold good;^X and Y direction fields. 

For jth column repeated section shown in Figure 4a, of the 
orthogonal cable network, the equations (2.20a) and (2.20b) 
developed in Section 2.1.3 can be rewritten as 



{FI 


3+t 


[T,J [T^] ' 

1 

Ur,} 1 
J 

. [T3] [T4] 


IlFl- 

L 3j 


(2.21a) 


or {ZL,, = [Tj i7J. (2.21b) 

J"** I J 

where [Tl is of order 2m x 2m, {"Hi {hi - and CF}.,., {Fl are 

3+1 j j+1 » j 

displacement and force vectors at right hand and left hand sides 
of the jth column respectively. 


2. 2. 1.1 TRANSFER MATRIX FOR LEFT END BOUNDARY OF CABLE NETWORK 

For element of cable network shown in Figure 4b, TM 
equation (2.4a) can be rewritten as 


1 . e. 


( h ' 



'N 

1 

(N 

< 


1 , 1 


*3 



= 

i 

X J 

.L 

*2^S, , 

i,L 

^i,1 

= 

A^T)„ + 

^i,1 

Vs 

^ ^i,L 


r) 


S. 


i,L 


(2.22) 


for i = 1 , 2, . . , ,m 


(2.23) 
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Fig. 4(a) column of orthogonal cable 




^m, L * ^m, I 

Fig. 4(b) Left end boundary , 
of the orthogonal 
network. 



Fig. 4(c) Right end boundary 
of the orthogonal 
network. 
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or. 


21i = [1^] IZIl 


(2.24) 


where and ^ Z are state vectors consisting of displace- 

ments and forces for all stations 3. , and S. ,, i = 1,2,...,m 

X ^ I ^3^ ^ 


and 


[Tl] = 


0 


0 

A. 


0 0 

0 0 

^3 ^ 

0 A, 

• • 

• • 

9 • 

0 0 

0 0 


0 

0 

♦ 

0 


0 

0 

* 

0 

A 


1 


0 0 
0 0 


A 3 0 

0 A, 


A^ 0 . 




0 0 
0 0 


A 4 0 

0 A. 


0 0 
0 0 


0 0 

0 0 

♦ • 
• • 
• • 

Ag 0 

0 A 

0 0 
0 0 


2 


A 4 0 

0 A 


4 


(2.25) 


2. 2. 1.2 TRANSFER MATRIX FOR RIGHT END BOWDARY OF CABLE NETWORK 

For element of cable network shown in Figure 4c, TM equa- 
tion (2.4a) can be rewritten as 


h‘ 

[ 

■ 81 

CM 

CQ 


' nA 

F 

J 

s. „ 

i,R I 

®3 

84 _ 


F . 

> 


^i,n+1 


i ♦ 




'i, R 

B R 

S 

^ i,n+1 

4- 

2 S. 
1 

i,R 

B.T) 

• ^i,n+1 

+ B.Fe 


i,n+1 

i,n+1 


(2.26) 


yfOP 1 ISS jIB 


* • • 


(2.27) 
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or. 


iZ}, 




(2.28) 


where iZl^^ and are state vectors consisting of displace- 
ments and forces for all stations and 3^ , i = 1,2, ...,m 

and 


LT^J = 


?1 0 

0 B, 

• • 

• • 

0 0 

0 0 


"3 " 

0 B, 


0 0 . 

0 0 . 


0 0 

0 0 

• ♦ 

« • 

• • 

0 

0 B 

0 0 

0 0 

• • 

• • 

• • 

B3 0 
0 B, 


1 


B 2 0 
0 B, 


0 0 
0 0 

B4 0 

0 B^ 

• • 

• • 

« • 

0 0 
0 0 


0 0 
0 0 


B^ 

0 

0 

0 


0 

B. 

j: 

0 

0 


®4 ^ 

0 B. 


(2.29) 


2. 2. 1.3 TRANSFER MATRIX FOR THE SYSTBl 


As mentioned earlier, the systsn is modelled to be 
combination of left end boundary, n repeated sections and right 
end boundary section. For each of these, TM has already been 
obtained. Combining these in accordance with the geometry of 
the system, one gets 

{Z}^ = [U] {Z}^ 

where [u] = [Tj^][tJ . . . [tJ[Tj^] 

= [Tr][T]" [Tl] 


(2.30) 
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Eauation (2.30) is the combined Tf'' relation between the state 
vectors on the right end and left end boundary of the cable 
network . 

2.2.2 TRANSFER MATRIX FOR NON-ORTHOGONAL CABLE NETWORK 

In this section the equations (2.21a) and (2.21b), are 

used to develop the TF'i for the non-ortho gonal cable network. Two 

types of non-orthogonal netw'orks are considered: 

. , i 

(i) A network in which the cables jtylnq in the Y direction 
start^ and end on X edges is referred to as type A non- 
orthogonal network, Figure 5. 

(ii) A network in which the cables lying in the Y direction 
start from X edge and end on Y edge and vice-versa, is 
referred to as type B non-orthogonal network, Figure 6. 

2. 2. 2.1 TYPE A NON-ORTHO(DNAL CABLE NETWORK 

Figure 5’ shows a square non-orthogonal type A cable 
network. The angle between the cables in the two directions is 

— 1 Le.njfti'- 

tan L/P, where L is the^-a-seae of the square network and P is the 
pitch between the cables. • In type A square network if there 
are m cables in the X direction, then there will be m+1 cables 
in the Y direction. For the cables lying in X direction, the 
portion of the cables between two successive joints and at 
boundaries are divided into two fields of equal length, for the 
cables lying in Y direction, the portion of the cables between 
two successive joints and at boundaries are treated as a single 
field. 

Equations (2.26a) and (2.26b), derived for the ortho- 

I 

gonal cable networks are applicable to the type A network. For 
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a set of cables to the left of a-a. Figure 5, let [l^i be the 
T''. The matrix [T^j is similar to the matrix [t] of equation 
(2.21b) with the difference tnat the stiffness matrix [k.^] used 
for the field in the X direction (left side of the first cable 
in the Y direction (j = 1)) is 



C O spX 
^ sinpX 


si xpX 

SiW 




sinpX 


cospX 

sinpX 


(2.31) 


where X — (m + 1 — i )/(ro + 1 ) 


(2.32) 


and i is for the ith cable in the X direction. Similarly CT2^ 
is the TM for cables to the right of section b-b. Figure 5. The 
stiffness matrix [K^j used for the field in the X direction 
(right side of the cable in the Y direction (j = m + 1 ))is 



cospX _ p 

P sinpX “ sinpX 

_ P o COS£X 

SinpX ^ sinpX 


(2.33) 


where X = i/(m + 1 ) 


(2.34) 


Keeping in view the geometry of the network the combined TM is 
obtained as 


= [t^JEt] . . . [T][T^]Cz}j. 

= Eu]Cz?j^ 

where Eu3 = Et^IEtI®"'' Et^I 


(2.35) 
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2. 2. 2. 2 TYPE 3 HON-ORTIDGOisiAL CABLE NET s/ORK 

In this chapter transfer matrix for the orthogonal and 
non-ortho gonal type A cable networks has been developed relating 
tno state vectors at left end and right end boundaries.. In case 
of type B non-orthogonal cable networks, the system is modelled 
in such a way that the combined TM relate the state vector of 
all the stations where Y-directed cable terminate and top and 
right hand boundaries to the state vector of all the stations 
from where Y-directed cable starts at bottom and left end bound- 
aries, Figure 6 , This is done in the following manner: 

Firstly state vector section S. as shown in 

bi 1 

Figure 7c, is obtained in terms of state vector ?Z)g at bottom 
boundary. The state vector IZ'^g at section can be obtained 
from state vector CZ?e by a TM which is generated by adding 
elements of TM for the Y-directed cable element between stations 
and $2 to that of a repeated section TM [t^] marked by 
thick lines. To proceed further, the state vector CzK ^t 

53 

section S 3 can be derived by substracting the terms for Y- 
directed cable element lying between the stations $2 
$2 g from state vector and then multiplying by TM [T^J and 

finally adding terms for Y-directed cable element as done 
earlier. Proceeding in the similar fashion TM relating state 
vectors for stations along top and right end boundaries, and 

left and bottom end boundaries can be obtained. This requires 

s 

TM for fields below and above any horizontal cable i which is 
obtainable from equation ( 2 . 6 a) as 


f 

n 


Cl 

^2 


Tl^ 

F 

d . 1 , 

C-:.. V 




' , F 


(2.36) 



24 



Fig. 7(a) Top boundary and last row of type-B nonorthogonal 
cable network. 



Fig. 7(b) i— and (i+1)— row of type-B nonorthogonal cable 
network. 



Fig. 7(c) Bottom boundary and first row of type B non- 
orthogonal cable network. 
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and 



respectively. It nay be noted that suffix J stands for a joint. 

The state vector can be expressed in terms of 

1 

state vector ?Z!g in the following way: 

IZ}^^ = [Tg] CZ^g (2.38) 

where [Tg] is TM for bottom end boundary. The TM ^Tg] is the 
same as expression equation (2.25) except that its elements are 
the elements of Tf/i of equation (2.36). Proceeding as stated 
above 


- [1122^ ^Zlg^ 


(2.39) 


where £Z}c: is state vector at 5^ 

^2 2 

Czig^ is state vector comprising of displacements and 
forces at stations g and bottom end boundary 



[TI^] = [T^][Tg] 


(2.40) 


(2.41 ) 


and 



(2.42) 


iri^j 


TI^(11) j TI^(12) 

“tT^TcT)-!” TI“(527 “ 


To arrive at state vector Czlc section two rows 

S3 3 

namely (m+1 ) and (2m+1 ) are eliminated from equation (2.39) in 
order to drop the Y-directed cable element between stations 


{Z!^^ = [1132) CZlg^ (2.43) 

It may be noted that [1132] can be obtained from [11^2^ equa- 
tion (2.40) by dropping the (m+1 )th row and last row. The 
order of the TH [TI32I is (2m) x (2m+2). Proceeding on the 
same lines, one arrives at 


^^^^m+1 ” ^^^2(m+1 )^ (2.44) 

where iZlc is state vector at section 

Vi 

fZ^LB state vector at stations comprising of 

bottom and left boundary except for station S . • 

m, L 

The order of Bi [i'i2(m+1 )'^ ^ (4m-2). 

To arrive at {Z}j for top boundary of the cable network, 

CZJg is multiplied by TM [t^] for the top boundary, thus 
m+1 

iZlr = CZL (2.45) 

^ ^ Vi 


= [T^] 


(2.46) 


or CZ}^ = [U’] {ZI^B 


(2.47) 


Denoting ?Z!jg as state vector for the set of stations 
at top and right hand boundary except ^ and for the 
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set of stations at left and bottom of the network except 
^n+1 , L 

= LUj (2.48) 

where [Uj is generated by adding elements of TK* of the right 
hand boundary fields (except for station to [U’ j . These 

are the elements which were dropped earlier from matrices 

i = 2,3, ...,m. 

2.3 CI-iARATT ERISTIC EQUATION 

The boundary conditions at left and right ends of the 
cable network are: 

= {0! and = (2.49) 


respectively. Substitution of these boundary conditions in the 
combined Tfv'i relation provides the characteristic equation. The 
combined TM relation for orthogonal and non-orthogonal type A 
cable network is 

CZl^ = [U] iZ}^ 


On partitioning, 


the displacements and forces^ om- 




i. 


I 

I 




(2.50) 


where [U^ ], [U2], [U3] and [U^] are square matrices of order m. 
Applying the boundary conditions of equation (2.49), one gets 



R 


u. 


1 I 


U. 


1 



( 2.51 ) 
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or [U 2 J {?'*L = 10 1 


(2.52) 


The natural frequencies of the cable networks are 
determined from zeros of the above characteristic equation. 
Similarly, the combined Tfi relation for non-ortho gonal type B 
cable network is 

On partitioning, one has 


( ^ ' 


— 

k = 

r J 

TR 


h I “2 


U, 


1 

ri 

F 

J 

>. - 


(2.53) 


LB 


^3 j '"4 

liere the boundary conditions at left and right ends, and top 
and bottom ends are 

= col, {■nl^=C0l, Ct) 1^ = {01 and Cnl^ =Col (2.54) 

respectively. Applying boundary conditions of equations (2.54), 
one gets 


1 ---' 



‘h 

r 

04 

D 

n 

” N 

0 

F 

. ■' 

t TR 


'«3 

' u, 

1 


F 

M J 


LB 


or 


[U 2 I 5F)j^b = 


( 2 . 55 ) 

(2.56) 


For its non-trivial solution, one must have 


IU 2 ' = ° 


the zeros of which lead to the natural frequencies. 
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CrLfi^PTER 3 

VIB' W LO N ANAL YS IS OF BEAi\ 'l GRILLAGE 


In the last chapter transfer matrices have been deve- 
loped for the transverse vibration of various types of cable 
networks. In this chapter FETM technique will be developed for 
the analysis of a beam grillage. As discussed in the case of 
cable network, orthogonal beam grillage also has the repeated 
characteristic. Figure 10a shows the repeated section of a 
orthogonal beam grillage. If the TM relating the state vector 
on left hand side and right hand side is known for a repeated 
section, the combined TM can be obtained, according to the 
geometry of the grillage. 

The grillage is analysed for its transverse vibration 
for (a) combined bending- torsional mode and (b) bending mode 
only by suitably idealizing the joints between the crossing 
member or adopting suitable stiffness properties. Rotary iner- 
tia and shear effects in the transverse vibration of the gril- 
lage are neglected,! . e. a beam element is idealized as an 
Euler* s beam. 

Firstly the TM® for the transverse and torsional vibra- 
tion of an individual (X- or Y-directed) member are written and 
then these are combined to from the TM of the repeated section 
of beam grillage. By multiplying the transfer matrices of 
repeated section and that of left end and right end boundaries, 
the combined TM [ul is generated. Using the boundary conditions 
the characteristic equation is arrived at which is solved for 


its roots 
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3.1 TRAHSFER MATRICES rOR INDIVIDUAL BEATv' ELH/ENT 

The governing equations of motion for transverse and 
torsional vibration of a beam element, shown in Figure 8a 
are [23] 


HI — 1 = 

m ’x" 

(3.1 ) 

3X^ 

at*^ 


and = 

p I 

P 3t^ 

(3.2) 

respectively, where 



E is modulus 

of elasticity 



I is area moment of inertia 

m is mass per unit length 

P is density 

C = 

G is shear modulus 


'^Ipc ~ torsional rigidity 


I is polar moment of inertia of equivalent circular 
pc 

section 

I is polar moment of inertia of rectangular section. 
Its TM for the transverse vibration is given by [l9] 


-W 




M 


V 


R 


LC. 


4S 

p r 


EI ^2 El ^3 


EI 1 EI ^2 


•^4 EI 4 H p P 

P l 2 ^2 P L ^3 ^o 

„4 EI « „4 EI n4 ^3 

P pC, P ^=2 P L 


LCj 


-W 


f 


M 


V 


(3.3) 
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v'here 3^ = 

(cosho + cosd)/2 

'■'1 “ 

(sinhp + sinp)/(2p) 

r* 

(coshp - cosp)/(2p^) 

-^3 ” 

(sinhp - sinp)/(2p ) 

4 

P 

mm^ lV(HI) 


(3.4) 


and is angular frequency of vibration. 

The Tl" for the torsional vibrations of the beam 
element, shovm in Figure 8a, is given by [l9] 

sinP 






T 

1 

I 


■R 


cosp 

L 

QIpc 

GI 

pc sinp 
“ L p 

cosP 



A 

e 


iJ 

J] 


(3.5) 


where 


|3 = cjL((I P)/(GI„^))''/2 


The grillage has beam elements in X and Y direction and as such 
their respective variables, material and geometric properties 
are identified by the suffixes X and Y. The equations (3.1) to 
(3.5) are non-dimensionalised with respect to the variables and 
material properties of the member lying in Y direction and are 
denoted by bar C) and are given below. Thus TM for Y direc- 
tion and X direction members take the shape as follows. 



Fig. 8(a) A beam element 



J j+i 

Fig. 8(b) X-directed beam element 



Fig. 8(c) Y- directed bean 
‘ element. 




12 3 j n-1 n 

Fig. 9 A mxn orthogonal beam grillage. 
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Transfer r/iatrix for Transverse Vibration of Y-directed Beam 
Slement: 


-W 

y 



“Co 

r 

^2 

^3 


-W 

y 

"y 



P '^3 

^0 

=1 

''^2 


¥ 

Y 

"y 



4 

P C2 

4 

n ^ 
p ‘^3 

=0 



1 “y' 


i +1 


p^'c, 


P^Cg 

^0 


<1 


where 'Vy = yL^, = (M^L^j/CE^l^) 

Vy = (VyL^VCEyly), = ( m ^ ( E ^ l ^ ) O.?) 

Transfer f/iatrix for Torsional Vibrations of Y-Directed Beam 


El ement: 





i+1 


r 

Vi '“''’I 


■\ 

cosp^ 


5 

3 

-r^p^sinp^ 

cosp^ 


w 


(3.8) 


where 


0 = 0 , T = (T L )/(E I ), 

y y* y '' y y'' ^ y y * 

p, = 0)((P I L^)/(G I 
^1 y py y'' y pcy 


r, = (I G )/{E I ) 
1 pcy Y^ Y Y 


(3.9) 
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Transfer j^iatrix for Transverse Vibrations of X-directed Beam 
Element: 



(3.10) 

where = (coshp2 + cosp2)/2 

= (sinhp2 + sinp2)/(2p2) 

D 2 = (coshp2 C0SP2)/(2P2) 

= (sinhp2 - sinp2)/(2p2) (3.11) 


w, = \/t^, t, = L^/L^, = (M^L^)/(E^I^) 

4 = 


Transfer Matrix for Torsional 
Element: 



Vibrations of X-directed Beam 



(3.13) 
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9^ = 9^, = (T^i^)/(E^I^) 

^3 = (^ocx°x>/‘Vy^ 


3.2 TR'V'JSFER ^'ATRIX FOR A REPEATED SECTION OF BEAfiS GRILLAGE IN 
COr.^.GINED BENDING-TORSIONAL VIBRATIONS 

The transfer matrix for a repeated section of a cable 
network has been developed in Section 2.1.3. On the similar 

To 

lines, the TIT for beam grillage can be developed, tor do so, one 
requires to deteimiine the stiffness matrix for a beam element. 

Using equations (3.6) and (3.8), state vector at 

section i+l , Figure 8c, can be written as in terms of state 
vector at section i as 



To express forces in terms of displacements, ecfuation (3.15) 


can be rewritten as 
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c? ^ 

y 


x+1 


"^121 I ^122 


t-A 




(3.16) 


Jx 


where 



fr 


' y ' 



’ and CF} = < 


— * 

-iv 



{f 


u tJ 


ID! 


Expanding equation (3.16) in terms of displacement and force 

vectors, force vectors IF ) and IF ! can be expressed in 

^ i+1 ^ i 

terms of displacement vectors ID J and ID ! as 


i+1 


y i 




1 

= 



^ ^ i+1 L 


[T 112 ] 


Ft i-i 
‘-^12'* 


or 


IF!, 


[^ 121 ] -['ri22^^^’^112^ C^i 22^ 1^^112^ 


ID!, 



fCD ) 1 


y i 




5Dv> 


^ l+l) 


(3.17a) 

(3.17b) 


Similarly, for a beam element lying in the X direction, equation 
( 3 . 17 ) can be rewritten as 


r-’4 = 

1 — 1 . 

04 

CM 

\ 


-CT 212 ]-' 

i 

’^°x^ 1 

1 j 

^ j+1 

^^221 

[7222^^212^ ''[T 211 ] 

^^222^^^212^ 

■ 

\ 

y 

Id ! 
"i+i'* 


(3.18a) 


or IF5^ = [K^e)] (3.18b) 

and are stiffness matrices of order (6x6) for X- 


and Y~directed beam elements respectively* 
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Figure 10a shows a repeated section of the beam grillage 
as for cable network and the numbering system adopted for fields 
in the section. The joint (i+m+l ) is shown separately to expl- 
ain the assemblage of i, i+m along the X direction and (2m+i), 
(2m+i+l) along the Y direction. For these fields the state 


vectors consisting of displacement vectors 






have been 


labelled with superscript of the number of element and subs- 
cripts of its two ends and are 
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Since for the repeated section having a number of fields and 
joints, the state vector consisting of displacements has to 
have a global numbering systems. In view of this above ment- 
ioned state vectors for elements meeting at (i+m+l), are desig- 
nated respectively as 
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m+1 



(a) A repeated section 



(b) (i + m+1)~ joint (c) Elements meeting at 

(i + m+1)t!^ joint 

Fig. 10 Finite element idealisation of a repeated section 
of beam grillage. 
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Having obtained the stiffness matrices for X- and 
Y-directed beam elements, these are assembled on lines similar 
to that of F©'’ as shown in equation (3.19). 


Qm+2 2m+2 ,2m+3 2m+4 3m+2 2^43 'in4 ^4^ ‘"»4 3fr42 



'’x2m43 

”x2nj44 
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Houation relates the forces and disolacement vectors in 

the finite element mesh of the repeated section of the beam 
grillage column as shown in Figure 10a. To transform equation 
(3.19) in the standard form of W, boundary conditions of nodes 
(m+1 ) and (?m+2) 

are aoolied. For the clamped case these are 




(3.20) 


where k =(m-(1)and (2m+2), 


Then equation (3.19) is 

reduced to the standard form of TMfA as done for the cable 
network i.e. 


»'3 


■ [T,] 

[Tgl ■ 


'"',1 


^ - 

_ P 3 I 

[’■4! 

- 

I '"', 


or IZ}^ = [T] 


(3.21a) 

(3.21b) 


where {DT and {F^ are displacement and force vectors respectively 
and 


[g] = [Tj] ([A^] - [AglrAjJ-l [A^]) 

[TjJ = [[AjJ [A3J-’ LA4JJ-' 

[T 3 ] = ([A 5 I - [A^][A 3 r' [a^]) [T,] - [A^IEAj]-' [Ag] 
CT4J = ([Ag] - [A^JLAj]"'' [A4J) [Tgl (3.22) 


Suffix 1 is used for left hand side of grillage finite element 
mesh and 3 is used for right hand side. 
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TT'(\?TSF^R ^'^TRTX for BRW." QRIILAGE WITH CLA.V 1 PED ENDS 

The TM for the system can be qenerated corresoonding to 
the geometry of the grid on lines similar to the cable network 
oroceeding from the left hand end, where 

= [T^JCz^L ^3.23) 

the state vector for (j+1 ) repeated section is given by 

= [Tjj [T^J (3.24) 

where j = 1 , . . . ,n. 

and finally for the right hand end as 
IZ}^ = [Tj^][Tf [T^J 

or IZ)^ = [U] CZ!^ (3.25) 

where [u] = [Tj^][Tf [Tj.] 

Relation (3.25) is the combined TM relation between the state 
vector on right end and left end of beam grillage. Note that 
the clamped conditions of the lower and upper horizontal ends 
have already been incorporated in the formulation. 

3.4 CHARACTERISTIC EQUATION 

For the beam grillage with clamped ends the character- 
istic equation can be derived using boundary conditions at left 
and right ends, similar to that of cable network. The boundary 
conditions for these ends are 

C = ^OJ and ~ 


(3.26) 
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v'fhore 
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oartitioning equation (3.?5), one has 
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V. J 


flow apDlying the above mentioned boundary conditions 


LU2J = iO) 


(3.27) 


(3.28) 


For its non-trivial solution one must have 
lU^l = 0 


(3.29) 


and the natural frequencies of the system are determined from 
the zeros of this equation. 


3.5 MODE SHAPES 

In Sections 3.2 and 3.3 the TM for the repeated section, 
left hand boundary and right hand boundary have been developed. 
These matrices will be used to develop the mode shape. 

Firstly, the homogeneous equation (3.28) is solved to 
compute Since is a null vector, so is completely 

known. Multiplying ^Zlj^ by [Tj^ 3, it gives i.e. state vectoi 

on the left hand side of 1st column. From ?Zl^, the displacement 
vector {D}^ can be found. Again by multiplying [t] (TM for 
repeated section), ^2^2 can be computed. In this way, the state 
vectors on the left hand and right hand side of the beam column 
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lying in the Y direction can be known. 

To find out the displacement vector CdIj at joints, 
at joints is Determined as following: 

Let be displacement vector at the left hand of the 

.ith column and the right hand of the jth column. Then 

by using eauation (3.19) 

+ [A4J (3.30) 

Knowing the displacement vector at the different sections and 
joints the mode shape can be plotted. 

In the case of repeated roots, the rank of equation 
(3.28) is reduced by two. This is solved by assigning values 
to two elements of the vector is quite likely that 

though two modal vectors are orthogonal to other modal vectors 
corresponding to the distinct roots, but are not mutually ortho- 
gonal. These are orthogonal! sed by standard procedure. 

Transverse vibration of beam grillage in pure bending 
is a special case of combined bending- torsion. As such the 
same procedure is adopted to find out the nodal frequencies and 
the mode shape by deleting terms corresponding to the torsional 
vibration. 
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CHAPTER 4 
RESUT.TS 

Tn Chapters 2 and 3, the TM® for the orthogonal and 
non-orthoqonal cable networks and orthogonal beam grillage 
have been developed. In this chapter, the method of solution 
IS presented so as to compute its natural frequencies. Modes 
for the beam grillage have also been found. The natural 
freouencies for various cable networks are compared with those 
obtained by TMM l21j and FEM l 7J , whereas those of the ortho- 
gonal beam grillage are compared with the ¥BA [18] and the 
plate analogy of beam grillage [18J. 

4.1 METHOD OF SOLUTION 

The natural frequencies of the system are determined 
from zeros of the characteristic equation lU^I =0* Newton- 
Raphson modified method is used to compute the zeros of charac- 
teristic equation. 

Pramod Gupta [21 ] used Golden Section Method for 
computation of natural frequencies. By using this method, he 
minimized the function. The Golden Section method is not 
fast converging and a large number of trials are conducted. 

For the fast convergence, Newton-Raphson modified 
method has been used in the present case. Since the function 
is not explicitly expressable, the Newton-Raphson method without 
modification cannot be used. To find out the tangent at a 
point consider two points, in its neighbourhood on its sides 
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nil'-' ^'l« lino joinim t\/o can be approximated as a tangent 

to f.ho polynomial function. It was conjectured that for a 
qivon error of convergence, this method does not miss the 

douM p T-ootc. 

loinputer orogrammes are developed for the computation 
of tie natural freguencies of networks and grillage and for 
irone saaoes for the grillage. 


4.:' SPnCIrllATION OF CABLE NETWORKS 


The Following data is used to calculate the natural 


Freouencies of the cable network: 

1 . Arm of the square cable network 

Cross-sectional area of the cable 

F. "/eight density of the material of 
cab] p 

4. Tension in cables lying in 
K direction 

5. Tension in cables lying in 
7 direction 


3.048 ra 
322.58 X 10“^ 
0.7855126 X 10^ N/m^ 
k 

333.673 XN 
k 

333.673 m 


4.3 RESULTS FOR ORTHOGONAL CABLE NETWORKS 

First four natural frequencies, important from the point 
of their design, of the orthogonal square networks (2x2 to 6x6), 
Figure 11, are calculated for the data given in Section 4.2 and 
are given in Table 1 . These values are compared with those 

obtained by F0i [7] and TMM [21 3. 

The data shows the closeness of the results. It is 

i 

worth pointing out that the present results are from one exact 
formulation and within the limit of computational errors. 
Whereas in the case of FEM-, the assumed displacement pattern 



(c) 4 X 4 (d) 5 X 5 (e) 6 X 6 

FIG 1 1 DTFFERLNT SQUARE ORTHOGONAL 

CABLE NETWORKS OF THE SANE 
SIDE LENbTH 
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Table 1 


I ‘ t lO f tonnon 
1 net' ,’o rl'. 


eG ^cpc) for the orthogonal 


sauare 


1 

’ a 

*" »/ ^ 

Fnt 


b. 

3x3 Net 


!(t)r'>nont 

T 

1 7 rr* ^ 

i oj 

I 

[ T?;r 

i [^‘j 

_i_ 

Node 

i FETM 

{(present 

{method) 

t 

i FSI 

I M 
« 

i B1M 

I 

'^1 1 

t i)r^ 

m 

Ml .6? 

58.91 

k 

1 1 

58.92 

60.45 

58.91 

1 . 


^7.43 

83.35 

^12 

90.74 

96.29 

90.70 



07,43 

88 . 35 

A 21 

90.74 

96.29 

90.70 


117.8b 

137.79 

117.81 

■^22 

117.85 

132.81 

117.81 


c . 

l/ode 

4x4 

let 


d. 

5x5 

Net 


[TzVu 
{(present 
{method ) 

j — 

{ F31 

i M 

1 

{ T1\1M 

1 M 

1 

Mode 

{ FETM 

{(present 

{method) 

1 

i FEM 

1 M 

t 

{ TMM 

1 P'J 

f 

^11 

■JO. 9 2 

59.88 

58.91 

^1 1 

58.92 

59.59 

58.91 

A 

12 

91 .69 

95.34 

91 .64 

'^1 2 

92.18 

94.73 

92.12 

^21 

91 .69 

95.34 

91.64 

^21 

92.18 

94,73 

92.12 

A 

^22 

117.85 

- 

117.81 

^22 

117.85 

“ 

117.81 


e. 6 x 6 Net 


f f 


Mode 

{ FETM 

; FEM 

{ TMM 

{(present 

{method) 

f 

[7j 

1 piJ 

f , - 

4i 

58.92 

59.41 

58,91 

^12 

92.46 

94.33 

92.40 

^21 

92.46 

94.33 

92.40 

.^22 

117.85 

124.72 

117.81 
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■> n ! :r, to o'^nr ^‘'•tiiration of stiffn^^ss and hence freauencies. 

iho .’Tror in Tssurod displacement pattern becomes more 
, jnooncf'd in (.he hlqher modes the error in the estimation of 
rrof'.i<^nc' es increase's with Lhe hi qher modes. Therefore, the 
n-rpnpnt method is betseT suited and more useful for computing 
r.,o ^’T'eouenci es and the modes. Further incorporating the FEM 
with t'le n.'v ttip formulation has been simplified, and the 
results ore qujte comoarable. 

4.4 "OR UOU-ORTHOGONAL CABLE NETWORKS 

In this section results for both Type A and Type B for 
non-ortho jonol sauare cable networks are presented. Natural 
frequencies for 2x3, 3x4, 4x5 and 5x6 Type A networks, shown 
in Figure 12, are given in Table 2 and compared, wherever poss- 
ible with those obtained by FEM [7j and TMM [2l] . Figure 13 
shows 2x3, 3x5, 4x7 and 5x9 Type B networks and their natural 
frequencies are given in Table 3. The results are comparable 
with those of FEM and TMM. 

4.5 SPECIFICATION OF BEAM GRILLAGE 

The following data is used to calculate the natural 
frequencies and mode shapes of the beam grillage: 

1 • Length of beam in X and Y direction : 

2. Area of cross section of beams in X : 
and Y direction 

3. Moments of inertia of the beams : 

lying in X and Y direction 

4. Elastic constants for the material : 
of beams lying in X and Y direction 


1 .27 m 

1 .93548 X 1 0""^ m^ 
1.04057 X 10"'^ 

207. 0 x 10^ N/m^ 



(•> 2x3 


(b) 3x4 






lED//: 
1 1 1 i i 



tc) 4 * S 

IG 12 


((J) S X 6 


dipflrgnf typf-a_ 
CABLG networks OF 
SIDE LENGTH 


nonorthogonal 
the same 



(a ) 2x3 


(b) 3x5 



(c) 4 X 7 (d) 5 X 9 

13 DIFPERENF TYPE -B NONOPTHOGONA 
CABLE NETWORKS OF THE SAME 
SIDE LENGTH 


Natural freauencies (cps) for Tyce 1 non-ortno' 



53 


. ' I'.q '!r>n<^itY o" thp ma tonal of 

'ar 


7 186 X 1 0 ^ Kg/m^ 


. n.t'ir o" tho material of : 89.0 x 10^ N/m^ 

in, 

-6 4 

V. iol of inertia of beams : 1 .040579 x 10 in 

lyin'! in 1 and Y direction 

^ A A 

c\, I'ol ir i:.oni-'‘nt of inertia of equi- : 0.3288034 x 10 m 

v'ilr’nt circular section of beams 
lying in X ind Y direction 


U ^ :LTS rOX ORTiiOaONAL BEAM dRILLAGE v/ITH CLAMPED ENDS 

4. ;.1 I -E: .’DIN 3 VIBRATIONS ONLY 

First four natural frequencies of the orthogonal beam 
grillage 2x2, 3x3, 4x4, 5x5 and 10x10, Figure 14, are calculated 
for the data given in Section 4.5 and are given in Table 4. 

These are compared with those obtained by FBI [18J and pia4re 
plate analogy of beam grillage [I8j. 

Comparison of the results shows that in this case the 
results obtained by FETM method are quite close to those obta- 
ined by FB^ and plate analogy. As discussed in the case of 
cable networks, the results obtained by FETM are exact within 
the computational errors. The present method has an advantage 
over FBVi and the plate analogy because of its simplified 
formulation. 

4.0.2 COMBINED BENDING-TORSIONAL VIBRATIONS 

First four natural frequencies for the clamped ortho- 
gonal beam grillage 2x2, 3x3, 4x4, 5x5 and 10x10, Figure 14, 
are calculated for the data given in Section 4.5 and are given 

in Table 5. 




Table 4 

Natural frequencies (cps) for the orthogonal beam grillage with clamped ends in bending vibrations only 
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Shnrrna i18j in the plate analogy and in +k 

FBfi assumed 

the beams to be torsionless i.e. the torsional < 54 . 

stiffness of 

the beam elements was considered to be negligible 1 

^ j ‘■Ae and zero 

which leads to bending vibration only. When the torsional 

stiffness (G) of the elanent is considered, the n , 

♦ ^ yrillage becomes 

more rigid as compared to that with negligible torsiof^ 3 j_ stiff 

, , one ~ 

ness and hence expects a higher value for the natural 

freauencies. Comparison of the data given in Tahi , , r- 

“J-es 4 and 5 

validates this justification. 


4.7 RESIXTS FOR ORTHOGONAL BEAJVl GRILLAGE WITH SImiitv 
SUPPORTED ENDS 


In this section, results for bending vibrations of 
orthogonal beam grillage with simply supported ends are presen 

ted. The natural frequencies for 2x2, 3x3 and 4 yA 

^ are given in 

Table 6 and compared with those of FEM [18J and 

P-J-ate analogy 

[18J. The results obtained are quite close to FEM and plate 

analogy. As discussed in the last section the present results 

are from the exact formulation and within the lim-; 4 . ^ 

Aiuit; of compu- 
tational errors. 


4.8 MODE SHAPES 


The mode shapes for the above mentioned cases hj*'^ been 
computed and some of these are shown in Figure 15 -(.q 
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Table 6 

Natural freauencies (cps) for the orthogonal beam grillage with 
aimply supoorted ends in bending vibrations only 
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Fig. 15 First four mode shapes for 3x3 clamped beam 
grillage in bending vibrations only. 





Fig. 16 First four mode shapes for 3x3 clamped beam 

grillage in combined bending torsional vibrations. 
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CHAPTER 5 

CONCLUSIONS 


In the present work, vibration analysis of various 
cable networks and beam grillages have been carried out with 
the help of FETM method. 

The natural frequencies corresponding to first four 
modes have been obtained for cable networks, orthogonal and non- 
orthoqonal. Also the natural frequencies and mode shapes for 
beam grillage (a) in bending vibrations and (b) in combined 
bending- torsional vibrations have been obtained. Comparison of 
the results with those reported in literature shows that the 
results are comparable. Further, the formulation of the pro- 
blem is obtained from the exact solution of its elements. 

In the vibration analysis of cable networks by FETM 
method the size of the TM is equal to twice the number of cables 
lying in one direction for orthogonal cable networks and non- 
orthogonal cable networks Type A, whereas in the case of non- 
orthogonal cable networks Type B, the size is less than twice 
the number of cables lying in both the direction by two. But 
in FEM the matrix size is dependent on the number of elements 
being considered. Therefore the matrix size in the present 
method is far lower than the FEM but is equal to the matrix 
size in TMM. In the present method, the analysis is quite 
simple as compared to TMM. So, the simplicity of analysis is 
the advantage of present method over TMM. 
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In vibration analysis of beam grillage by the present 
mothod the size of TM is also small as compared to the FB^. 

Based on the present work it can be concluded that 
F5T?'5 method is quite powerful and economical tool for the vibra- 
tional analysis of systems which have repeated nature. 
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